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Radiative decays and the nature of heavy quarkonia 
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We argue that the photon spectra in radiative decays of various heavy quarkonium states provide 
important information on their nature. If two of these states are in the strong coupling regime, we are 
able to produce a parameter-free model independent formula, which holds at next-to-leading order 
and includes both direct and fragmentation contributions. When the formula is checked against 
recent CLEO data it favors T(2S) and T(3S I ) in the strong coupling regime and disfavors T(IS') in 
it. 
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Heavy quarkonium systems have played a major 
role in our understanding of QCD (see for a re- 
view). Inclusive heavy quarkonium decays to light 
particles involve a short distance process, in which 
the heavy quark and antiquark annihilate into glu- 
ons or photons, and a long distance process, which 
takes into account that the heavy quark and anti- 
quark are bound in a color singlet state. Due to the 
asymptotic freedom of QCD the short distance pro- 
cess can be calculated using perturbation theory in 
a s (m), m being the heavy quark mass. It was be- 
lieved for some time that in order to have a color 
singlet state the heavy quark and antiquark in it 
had to be in a color singlet state themselves. This 
allowed to parameterize the long distance process by 
a wave function at the origin (or derivatives of it). In 
spite of the fact that this wave function was not com- 
putable from perturbative QCD (pQCD), it canceled 
in suitable ratios, and hence predictions depending 
on pQCD only could be put forward. Nevertheless, it 
was soon noticed that the above framework had diffi- 
culties due to the infrared divergences which showed 
up at higher orders of a s (m) |2||. The introduction 
of color octet operators in the framework of Non- 
relativistic QCD (NRQCD) provided a rigorous 
way to understand, and to deal with, these diver- 
gences . The price to be paid was the introduction 
of new non-perturbative parameters, the color octet 
matrix elements, which made predictions depending 
on pQCD only more difficult to obtain. 

Later on, it was pointed out that of the hierar- 
chy of relevant scales in heavy quarkonium m 3> 
mv tow 2 , v being the relative velocity, only the 
first inequality was exploited in NRQCD and that 
the velocity counting proposed in |3| did not nec- 
essarily follow from it. It was also shown how the 
last inequality could be exploited by constructing a 
further effective theory, namely Potential NRQCD 
(pNRQCD) 0,0 (see la for a review). The inter- 
play of Aqcd with the scales mv and mv 2 above 



dictates the degrees of freedom of pNRQCD. Two 
regimes have been identified : (i) the weak coupling 
regime, Aqcd mv 2 , and (ii) the strong coupling 



regime mv <C Aqcd ;$ 



In the weak coupling 



regime, the relevant degrees of freedom are a singlet 
and an octet wave function fields interacting with 
(perturbative) potentials and with ultrasoft gluons. 
The original NRQCD counting belongs to this 
regime. Moreover, the use of diagrammatic and 
quantum mechanical methods allows to carry out 
explicit calculations at higher orders of a s (mv) 0. 
In the strong coupling regime, the essential degree of 
freedom is a singlet wave function field interacting 
with a (non-perturbative) potential. In this regime, 
the matrix elements of both color singlet and color 
octet NRQCD operators reduce to wave functions at 
the origin plus a few universal parameters . 

Due to the fact that none of the scales involved 
in the hierarchies above are directly accessible ex- 
perimentally, given a heavy quarkonium state, it is 
not obvious to which regime it must be assigned. 
Only the T(1S) appears to belong to the weak cou- 
pling regime, since weak coupling calculations in 
a s (mv) converge reasonably well. The fact that the 
spectrum of excitations is not Coulombic suggests 
that the higher excitations are not in the weak cou- 
pling regime, which can be understood from the fact 
that O(Aqcd) effects in this regime are proportional 
to a high power of the principal quantum number 
[lfl ]. Nevertheless, there have been claims in the 
literature, using renormalon-based approaches, that 
also T(2S') and even T(3S) can also be understood 
within the weak coupling regime [ill ]. We argue in 
this letter that the photon spectra in semi-inclusive 
radiative decays of heavy quarkonia to light hadrons 
provide important information which may eventu- 
ally settle this question. 

The contributions to the decay width of a state n 
can be split into direct and fragmentation, dT n /dz = 
dT^ r / dz+dT^' ag / dz. Direct contributions are those 
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in which the observed photon is emitted from the 
heavy quarks and fragmentation contributions those 
in which it is emitted from the decay products (light 
quarks), z 6 [0,1] is defined as z = 2E~ f /M n (M n 
is the mass of the heavy quarkonium state), namely 
the fraction of the maximum energy the photon may 
have in the heavy quarkonium rest frame. The ap- 
proximations required to calculate the direct contri- 
butions are different in the lower end-point region 
(z — > 0), in the central region (z ~ 0.5) and in the 
upper end-point region {z — ► 1) of the spectrum 0] • 
We shall restrict our discussion to z in the central re- 
gion, in which no further scale is introduced beyond 
those inherent of the non-relativistic system. Con- 
sequently, the photon spectrum can be expressed in 
terms of matrix elements of local NRQCD operators 
Q with matching coefficients C[Q](z) which depend 
on m and z. 



dz 
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(Q)r 



(1) 



Sq is an integer which follows from the dimension of 
Q, (Q) n := (V Q (nS)\Q\V Q (nS)) and V Q {nS) stands 
for a vector S*-wave state of principal quantum num- 
ber n. The fragmentation contributions read |l3| 



Q 



a=Q,Q,9 Q 
(Q)r 



(2) 



D ai (x, m) are the fragmentation functions and 
C a [Q](x) the partonic kernels. It is important for 
what follows that the /q(z) are universal and do 
not depend on the specific bound state n. Due to the 
behavior of the fragmentation functions above, the 
fragmentation contributions are expected to domi- 
nate the spectrum in the lower z region and to be 
negligible in the upper z one. In the central region, 
in which we will focus on, they can be treated as a 
perturbation, as we will show below. 

Let us first consider the weak coupling regime, for 
which the original NRQCD velocity counting holds 
0. The direct contributions are given at leading 
order (LO) by the 0\ ( 3 <Si) operator; the next-to- 
leading order (NLO) (v 2 suppressed) term is given 
by the V\ \S\) operator. The contributions of color 
octet operators start at order v A and are not a" 1 (to) 
enhanced in the central region. The fragmentation 
contributions are more difficult to organize since the 
importance of each term is not only fixed by the ve- 
locity counting alone but also involves the size of 
the fragmentation functions. It will be enough for 



us to restrict ourselves to the LO operators both 
in the singlet and octet sectors. The LO color sin- 
glet operator is 0\ ( 3 <Si) as well. The leading color 
octet contributions are v 4 suppressed but do have a 
a~ x (m) ~ 1/v 2 enhancement with respect to the sin- 
glet ones here. They involve C 8 ( 3 <Si), O s (^o) and 
8 ( 3 P ) • Then in the central region, the NRQCD 
expression (at the order described above) reads 



dz 
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<08( 3 Si)>. 



+ fo S ( 3 Pj){ Z ) 



fc 
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C 1 [ 3 S 1 ](z) and C[ [ 3 5i] (z) are the only short dis- 
tance matching coefficients that will be eventually 
needed. They can be found in 0] and respec- 
tively. If we are in the strong coupling regime and 
use the so called conservative counting, the color 
octet matrix elements are suppressed by v 2 rather 
than by v 4 . Hence we should include the color octet 
operators in the direct contributions as well. In prac- 
tise, this only amounts to the addition of Cs's to 
the / 0g 's. Furthermore, f 0l ( 3 Si) ( z ), /o 8 ( 1 So)( z ) and 
fo s ( 3 Pj)( z ) are proportional to D gl (a;, to), which is 
small (in the central region) according to the widely 
accepted model jj^- /o s ( 3 5i)( z ) is proportional to 
D ql (x,m), which has been measured at LEP [l^. It 
turns out that numerically fo^St) ( z ) ~ Cs[ 3 Si](z) 
in the central region. Therefore, all the LO frag- 
mentation contributions can be treated as a pertur- 
bation. Consequently, the ratio of decay widths of 
two states with different principal quantum numbers 
is given at NLO by 

dT n 

(0 1 ( 3 ^ 1 ))» c' 1 [*s 1 }{z)n™ {3Si) 

dT r (0 X (35l)) P V + Ci[ 3 Si](z) TO 2 



dz 



/o 8 ( 3 5 , l)( Z )^ nr | /Os( 1 g )( 2: )^ nr 

fo s ( 3 Pj)(z) n7 O s { 3 P ) 



C\ [35i] (z) m 2 



where 
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Note that the a s (m) corrections to the matching 
coefficients give rise to negligible next-to-next-to- 
leading order (NNLO) contributions in the ratios 
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above. No further simplifications can be achieved 
at NLO without explicit assumptions on the count- 
ing. If the two states n and r are in the weak cou- 
pling regime, then T^p 1 ^s 1 ) = fn(E n — E r ) (E n is the 
binding energy) [li}. In addition, the ratio of matrix 
elements in front of the rhs of J3J can be expressed 
in terms of the measured leptonic decay widths 



(0 1 ( 3 S 1 )) n _ T(V Q (nS) 



(Oi( 3 5i)) r T(V Q (rS)^e+e-) 



1- 



lmg ee ( 3 5i) E n - E r 
Im/ ee ( 3 5i) m 



(6) 



Im(; ee and Im/ ee are short distance matching coef- 
ficient which may be found in 3]. Eq.© and the 
expression for T^p^s-A a ^ so if both n and r are 
in the strong coupling coupling regime 0, 0] , but 
none of them does if one of the states is in the weak 
coupling regime and the other in the strong coupling 
regime. In the last case the NRQCD expression de- 
pends on five unknown parameters, which depend 
on n and r. If both n and r are in the strong cou- 
pling regime further simplifications occur. The ma- 
trix elements of the color octet NRQCD operators 
are proportional to the wave function at the ori- 
gin times universal (bound state independent) non- 
perturbative parameters [H, El- Since (0i( 3 5i))„ 
is also proportional to the wave function at the ori- 
gin, the latter cancels in the ratios involved in 
Hence, Wq = for the octet operators appearing in 
(@J. Then, the pNRQCD expression for the ratio of 
decay widths reads 



dz 
dT r 

dz 



(0i( 3 Si))» 

(Ol( 3 Sl))r 



X 1 + 



C[ [ 3 5i] (*) 1 
Ci [ 3 Si] (z) m 



(E n -E r )\ (7) 



Therefore, in the strong coupling regime we can pre- 
dict , using pNRQCD, the ratio of photon spectra at 
NLO (in the v 2 , (A QCD /m) 2 9] and a s (^/m A QC d) X 
y/A-QCT)/m 01 expansion), which is the main result 
of this letter. On the other hand, if one of the states 
n is in the weak coupling regime of pNRQCD, TZq 
will have a non-trivial dependence on the principal 
quantum number n and hence it is not expected to 
vanish. Therefore, expression J7J provides invalu- 
able help for identifying the nature of heavy quarko- 
nium states. If the two states are in the strong cou- 
pling regime, the ratio must follow the formula (|7|): 
on the other hand, if (at least) one of the states is 
in the weak coupling regime the ratio is expected 



to deviate from JjJ, and should follow the general 
formula (J2J. We illustrate the expected deviations 
in the plots (dashed curves) by assigning to the un- 
known IZs in (J3J the value v 4 (v 2 ~ 0.1), according 
to the original NRQCD velocity scaling. 

Recently CLEO [2(J has measured the photon 
spectrum for the Y(15) radiative decay (with very 
good precision) and also (for the first time) the spec- 
tra for the T(25) and Y(35) radiative decays. We 
will use this data to check our predictions. In order 
to do the comparison we use the following proce- 
dure. First we efficiency correct the data (using the 
efficiencies modeled by CLEO). Then we perform the 
ratios 15/25, 15/35 and 25/35 (we add the errors of 
the different spectra in quadrature) . Now we want to 
discern which of these ratios follow eq. J7J) and which 
ones deviate from it; to do that we fit eq.© to each 
of the ratios leaving only the overall normalization 
as a free parameter (the experimental normalization 
is unknown). The fits are done in the central re- 
gion, that is z e [0.4,0.7], where eq.Q holds. A 
good (bad) x 2 obtained from the fit will indicate 
that the ratio does (not) follow the shape dictated by 
eq.0. In figures HJ El andCHwe plot the ratios 15/25, 
15/35 and 25/35 (respectively) together with eq.Q 
and the estimate of Q mentioned above (overall 
normalizations fitted for all curves, the number of 
d.o.f. is then 45). The figures show the spectra for 
z £ [0.2, 1] for an easier visualization but remember 
that we are focusing in the central z region, denoted 
by the unshaded region in the plots. The theoreti- 
cal errors due to higher orders in a s (m) and in the 
expansions below (JJJ) are negligible with respect to 
the experimental ones. For the 15/25 ratio we ob- 
tain a X 2 /d-°-f-|iS/2,s ~ 1.2, which corresponds to an 
18% CL. The errors for the T(35) photon spectrum 
are considerably larger than those of the other two 
states, this causes the ratios involving the 35 state 
to be less conclusive than the other one. In any case 
we obtain x 2 /d.o.f.|is/3s ~ 0.9, which corresponds 
to a 68% CL, and x 2 /d-°-f-|2S/3S ~ 0.75 which cor- 
responds to an 89% CL. Hence, the data disfavors 
T(15) in the strong coupling regime but is consis- 
tent with T(25) and T(35) in it. For completeness 
we also quote the numbers corresponding to the es- 
timate of eq. 10} (dashed lines): x 2 /d.o.f.|is/2S ~ 
1.93 (.02% CL), x 2 / d -°- f -li5/35 ~ -56 (99% CL) and 
X 2 /d.o.f.| 2S/3S - .66 (96% CL). Recall that these 
curves are only estimates to illustrate what the dif- 
ferences from eq. Q to eq. (J2J may be and do not 
intend to best fit data. 

In summary, using pNRQCD we have worked 
out a model-independent formula which involves the 
photon spectra of two heavy quarkonium states and 
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FIG. 1: Ratio of the T(15) and T(2S) photon spec- 
tra. The points are obtained from the CLEO data |20)|. 
The solid line is eq.Q (overall normalization fitted), the 
dashed line is the estimate of £IJ (see text). Agreement 
between the solid curve and the points in the central (un- 
shaded) region would indicate that the two states are in 
the strong coupling regime. 



crease of the current experimental errors for T(2S) 
and, specially, for the T(35) is necessary to confirm 
this indication. This is important, since it would val- 
idate the use of the formulas in [l9| , and others 
which may be derived in the future under the same 
assumptions, not only for the T(2S) and T(3S I ) but 
also for the %b(2P)s, since their masses lie in be- 
tween, as well as for their pseudoscalar partners. 
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FIG. 2: Same as fig. 1 for T(1S) and T(3S) . 



holds at NLO in the strong coupling regime. When 
this formula is applied to the Upsilon system, cur- 
rent data indicate that the T(2S) and the Y(3S) are 
consistent as states in the strong coupling regime 
whereas the T(IS') in this regime is disfavor. A de- 
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FIG. 3: Same as fig. 1 for T(2S) and T(35) . 
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